Slip line growth as a critical phenomenon 
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We study the growth of shp hne in a plasticaUy deforming crystal by numerical simulation of a 
double-ended pile-up model with a dislocation source at one end, and an absorbing wall at the other 
end. In presence of defects, the pile-up undergoes a second order non-equilibrium phase transition 
as a function of stress, which can be characterized by finite size scaling. We obtain a complete set of 
critical exponents and scaling functions that describe the spatiotemporal dynamics of the slip line. 
Our findings allow to reinterpret earlier experiments on slip line kinematography as evidence of a 
dynamic critical phenomenon. 
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Plastic deformation in crystals is due to the motion 
of dislocations driven by the external stress. In re- 
cent years, experimental and theoretical work has shown 
that dislocation dynamics is a complex intermittent phe- 
nomenon involving the collective motion of many inter- 
acting dislocations 1]. In particular, deformation test on 
micron-scale crystals has revealed intriguing size effects 
and power law distributed strain bursts ^, iSj , which have 
been reproduced in dislocation dynamics simulations 
Intermittent dislocation motion is a generic feature of 
plasticity, not only in micron-scale samples, as shown re- 
cently by acoustic emission measurements 0] and by 
earlier reports from slip line kinematography in macro- 
scopic samples 0, iS,, ^] . These experimental observation 
lead to the idea that plastic yielding is a non-equilibrium 
critical point [ll], similar to the jamming transition ob- 
served in soft and glassy materials ,10J or the depinning 
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transition for disordered elastic manifolds 

The yielding transition has been investigated in vari- 
ous dislocation models, from the dynamics of an individ- 
ual flexible dislocation interacting with quenched random 
impurities such as solute atoms [l^ [l3| , to the dynamics 
of several rigid dislocations moving on single slip systems 
3, 1^ • While these models provided a good understand- 
ing of the yielding transition in simplified conditions, less 
is known about the role of dislocation nucleation and mul- 
tiplication for the critical behavior. In very clean single 
crystals, with a very small initial dislocation density, dis- 
locations are most likely nucleated from sources present 
at the surface of the sample, where is easy to find defects 
(steps, scratches) acting as stress concentrators. In this 
case, the onset of plasticity corresponds to the creation 
and propagation of slip lines through the entire cross sec- 
tion of the crystal. 

A slip line can be envisaged as a queue of dislocations, 
a pile-up, pushed through a series of obstacles, such as 
solute atoms or immobile dislocations from other glide 
planes (see Fig. [1]). Experimentally slip lines can ter- 



minate or propagate depending on the value of the shear 
stress, temperature, crystal structure and types of defects 
[l6| . A transition from homogeneous to inhomogeneous 
slip, with increasing impurity concentration, is observed 
experimentally in fee alloys 17, Here, we investi- 

gate the dynamics of a double-ended pile-up in presence 
of defects, with a source of dislocations at one end and 
an absorbing wall at the other. In the model, the stress 
dependence of the stationary dislocation density, velocity 
and strain rate obey finite-size scaling. This result indi- 
cates that the transition observed in earlier experiments 
should be reinterpreted as a signature of a second-order 
non-equilibrium critical phenomenon. The scaling expo- 
nents we measure in the present model are different from 
those found in the corresponding homogeneous system, 
where nucleation is not considered and the dislocation 
density is constant [19]. 

We consider the pile-up as a group of identical edge or 
screw straight dislocations parallel to the z axis that can 
move in the positive x direction in the plane j/ = when 
the net force acting on them is positive. This corresponds 
to an effective one-dimensional model in which the dis- 
locations are generated from a source in the left side of 
a line of length L (see Fig. [T]). They interact with each 
other and with a disordered stress landscape provided by 
solute atoms, or other defects, and disappear when they 
reach the right side of the line. The dislocations have co- 
ordinates Xi, with i = \, ...N , where N — N{t) depends 
on time t. The dislocation at i = 1 is immobile and rep- 
resents the source, as we discuss below. The dislocations 
for i = 2, ...N are mobile and have constant Burgers vec- 
tor bi = b. The Burgers vector is directed along x for 
edge dislocations and along z for screw ones. 

To describe the dynamics of mobile dislocations we use 
an over-damped equation, so that the velocity of dislo- 
cations depends linearly on the resolved shear stress ex- 
erted on it [iBi- The equation of motion for the mobile 
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FIG. 1: At the top, we show a three dimensional sketch of 
a shp plane line containing a pile-up of dislocations emitted 
by a source placed on the left surface and absorbed at the 
right hand side of the sample. Our simplified model is shown 
in the bottom part of the figure: dislocations move on a line 
of length L and interact with point impurities. The source 
is modeled as an immobile dislocation with varying Burgers 
vector placed at x\. In order to take into account the effect 
of open boundary conditions, interaction stresses include the 
contribution from an infinite series of image dislocations. 



dislocations is given by 
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where x is an effective viscosity and a is the external 
stress. The interaction stress cr-"* between dislocations i 
and j is computed taking into account the image stresses 
of dislocation j due to the open boundary conditions, 
while (7™^ is due to the interaction between dislocation 
i and its own images. A compact expression of the inter- 
action and image stresses can be obtained by performing 
the sum over the images [13 j yielding 
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where /i is the shear modulus, fc = tt for screw disloca- 
tions and k = 2ti{\ — v) for edge dislocations, v is the 
Poisson ratio. We notice here, that the sum over the im- 
ages is exact only in the case of screw dislocations. For 
edge dislocations, there is an additional subdominant cor- 
rection scaling as 1 /r^ that we neglect here since it should 
not influence the scaling behavior. The last term in Eq.[T] 
represents the interactions with pinning centers placed 
at randomly chosen positions Xp with P = l,...,Np. 



The detailed shape f{x) of the individual pinning force 
is inessential for most purposes, provided it is of short- 
range nature, and in this case it is given by 
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where ^p is the range of the interaction and /o controls 
its strength. 

Dislocations are typically generated by Frank-Read 
like sources, which can only be represented in a three 
dimensional model. In lower dimensions, it is custom- 
ary to model the source phenomenologically by creating 
dislocations with a certain rate. The drawback of this 
approach is that the new dislocation produces an arti- 
ficial discontinuity in the stress field. To overcome this 



problem, we employ a method suggested by Zaiser 21 1 
in which a source is represented by an immobile disloca- 
tion, placed at position xi, with a time-dependent Burg- 
ers vector bi{t) growing with stress. When bi{t) — b, a 
new mobile dislocation is emitted from the source whose 
Burgers vector is reset to zero. The evolution equation 
for bi (t) is given by 
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where 9 is the Heaviside step function and xi is a damp- 
ing constant that we set equal to xi = x/b- The effective 
stress al^^ = a + cr™* -|- cr™^ is the sum of the constant 
external stress tr, the stress cr™* produced by the inter- 
action between the source and the mobile dislocations 
(including the relative images), and the stress cr^™^ pro- 
duced by the interaction between the source and its own 
images. These stresses are obtained from Eq. [2] observing 
that ar* = Ef=2 and - cr^f . 

Integrating numerically Eqs. [T] and 01 we analyze the 
dynamics of the pile-up as a function on the external 
stress a and the system size L. The units of time, space, 
and forces are chosen so that & = 1, x = 1 and /i/fc = 1. 
For the simulations reported here, we considered parame- 
ters L = 256, 512, 1024, 2048, 4096, 8192, xi = 16 and the 
pinning centers are Poisson distributed with an average 
spacing dp = L/Np — 2 with /o = 1 and ^ = 1. From the 
experimental point of view, a key quantity describing the 
growth of the slip line is the plastic strain rate e given by 
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where p ~ N/L is the dislocation density and v — 
J2 ii is the average dislocation velocity. Notice 
that all these quantities are defined per unit dislocation 
length, given the effective one dimensional geometry of 
our model. 

Since the strain rate is simply the product of the dis- 
location density and average velocity, we study directly 
these two quantities. We find that after an initial tran- 
sient the density and the velocity reach a steady state 
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value {ps and Vs respectively) that depends on the sys- 
tem size L and the applied stress a as shown in the inset 
of Fig. [51 The graphs are suggestive of a non-equilibrium 
phase transition controlled by the stress between a pinned 
phase at low stress and a moving phase at large stresses. 
The curves become sharper close to the depinning point 
as L is increased, as expected when finite-size effects are 
present. To confirm this idea we perform a scaling col- 
lapse according to 
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where the scaling function f{u) fulfills the limits 



if u < 1, 
if u > 1, 



and for g{u) they are 



if M < 1, 
if M > 1. 
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The best collapse is obtained using Uc — 1.05 ± 0.05, 
V = 2.85 ± 0.05, a/v = 0.35 ± 0.02 and (3/v = 0.17 ± 
0.02, as shown in Fig. [2l These exponent combinations 
correspond to a = 1.00 ± 0.02, [i = 0.48 ± 0.02, and yield 
a scaling form for the strain rate of the type 
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where h{u) = f{u)g{u), as we have also verified directly. 

We have also analyzed the slip line growth dynamics 
in the transient regime. The time dependence of the 
dislocation density and velocity can also be characterized 
by finite size scaling functions which at the critical point 
(T = (Tc are given by 



pit,L) - L-^/-^ft[t/L% 
v{t,L) - L~f>/''gt[t/L% 
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where z the dynamic exponent. The best collapse is ob- 
tained for z = 1.25 ± 0.02 as shown in Fig. [3l The scal- 
ing functions, for small values of the argument, scale as 
ft{u) ~ u'^, with C = 0.55 ± 0.05 and gt{u) ~ u^^ , with 
6 — 0.10 ± 0.05. Hence, the strain rate in the initial 
phase grows as e{t) ~ Notice that the scaling ex- 

ponents are considerably different from what is expected 
for a regularly spaced pileup with periodic boundary con- 
ditions, where the density is constant (hence a = 0) and 
the critical exponents are (3 ~ 0.78, z ~ 0.78 and ly ~ 1.5 

To elucidate the role of the boundary condition and 
characterize the internal morphology of the pile-up, we 
report in FigH] the stationary density ps{x,L) and ve- 
locity Vs{x,L) profiles for different values of the system 
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FIG. 2: (a) The average stationary density of the pile-up as a 
function of the system size L for different values of the applied 
stress a (inset). The data collapse is consistent with the scal- 
ing hypothesis in Eq[6]with the exponents a = 1.00±0.02 and 
l/v — 0.35 ± 0.02. (b) The average stationary velocity of the 
pile-up as a function of the system size L for different values 
of the applied stress u (inset). The data collapse is consis- 
tent with the scaling hypothesis in Eq|S]with the exponents 
f3 = 0.48 ± 0.02 and u = 2.85 ± 0.05. 
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FIG. 3: The evolution of the dislocation density and velocity 
(inset) as a function of the system size L at the critical point 
a — ac — 1.05. The best collapse is obtained for a/u = 
0.35 ± 0.02, P/u = 0.17 ± 0.02 and z = 1.25 ± 0.02 which is 
consistent with the scaling collapse in Eq. (2] 
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FIG. 4: The density ps{x, L) and velocity Vs (x, L) (inset) pro- 
file for different values of the system size L, for a = 1.05 
(critical regime). In the bulk of the line we have the power- 
law behavior ps{x,L) ~ x^"' jL^ and Vs{x,L) ~ x'^/L'^ with 
7 = 0.25 ± 0.02, V = 0.09 ± 0.02 and (j> = 0.42 ± 0.02. These 
behaviors are consistent with the scaling of the stationary 
density and velocity in the critical regime. 



size L, for a — <7c — 1.05. We observe inhomogeneities for 
both density and velocity profiles which can be described 
as power laws: Psix,L) ~ x-'^/L'^ andvs{x,L) - x^/L"^, 
with 7 ~ 0.25±0.02, ^ ^ 0.09±0.02 and (j) = 0.42±0.02. 
Notice that the strain rate profile is approximately con- 
stant since the two power law cancel out in the product. 
This behavior is consistent with the scaling of the sta- 
tionary density and velocity as described in Eq. [S) In 
particular, we have that 



Ps{(Tc,L) = y[ Ps{x,L)dx^L' 

^ Jxi 
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Hence, we have the scaling relation 7 -I- '0 = ol/u, that is 
verified by the numerical values of the exponents. Simi- 
larly, from the steady-state velocity equation 

..(..,L)=^-'^"t'^^"-^"'^^'"^L-(^-^) (12) 

we obtain the relation (j) — j = (3/v, which is again in 
agreement with our numerical estimates. 

In conclusion, we have studied the slip line formation 
at the initial stage of plastic deformation in a crystal 
by means of the double-ended pile-up model finding that 
in presence of pinning centers (quenched disorder) the 
model exhibit a non-equilibrium phase transition. As 
a consequence of this dislocation density, velocity and 
strain-rate are described by finite size scaling. Finite size 
scaling has direct implications for size effects: the size de- 
pendence of the yield stress ay observed in micron scale 



plasticity Q . Considering the scaling law in Eq. [6l we 
expect that the yield stress for finite L grows towards the 
asymptotic value CTc according to (Ty(L) — — A/L^/" , 
where ^ is a positive constant. This type of inverse size 
effect, with the strength increasing with the sample size, 
is due to the larger back-stress exerted on the source by 
the pile-up as its length is increased. In more general 
cases, involving many sources and several slip lines, the 
constant A is expected to be negative as shown in other 
models of the yielding transition [ij . 
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